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Abstract: Within the frame of cosmologies where Dark Energy (DE) is a self-interacting 
scalar field, we allow for a CDM-DE coupling and non-zero neutrino masses, simulta- 
neously. In their 0-0 version, i.e. in the absence of coupling and neutrino mass, these 
cosmologies provide an excellent fit to WMAP, SNIa and deep galaxy sample spectra, at 
least as good as ACDM. When the new degrees of freedom are open, we find that CDM-DE 
coupling and significant neutrino masses 0.1 eV per u species) are at least as likely as 
the 0-0 option and, in some cases, even statistically favoured. Results are obtained by 
using a Monte Carlo Markov Chain approach. 
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1. Introduction 

In the late 90's Hubble diagrams of SNIa [||] became sufficiently precise to allow the unex- 
pected conclusion that the cosmic expansion is accelerated. This agreed with fresh Cosmic 
Microwave Background (CMB) Q and large scale structure (LSS) ^ data, strongly in- 
dicating that the background metric is spatially flat, while the matter density parameter 
Qo^rn — 0.27 is much below unity. It was then natural to infer that the rest of the energy 
budget (up to Qo — 1) is responsible for the cosmic acceleration; dubbed Dark Energy 
(DE), it ought to be a smooth non-clustering component, with a state parameter close to 
-1. 

All data available up to now can be accomodated in a cosmology where DE has a state 
parameter w = —1, a model equivalent to introducing Einstein's cosmological constant. 
This minimal model is dubbed ACDM (or cosmic concordance) cosmology. However, DE 
nature has yet to be properly understood. 

A component with w = —1 could be false vacuum. If vacuum energy does not vanish, 
its expected density is [mp-. the Planck scale), and the measured density implies a fine 
tuning ~ 1 : 10^'^^. But, even referring to the last phase transition, supposed to occur when 
the cosmic temperature was Tew ~ 100 GeV, the fine tuning is still ~ 1 : 10^^. 

The coincidence paradox is even more severe. If one does not want to indulge to an- 
thropic perspectives, a vacuum energy level just allowing structure formation, and stopping 
the process just when it has completed, can be hardly accepted without a justification. 

Alternatives to false vacuum were then proposed, aimed to avoid fine tuning and 
coincidence. They mostly bring new parameters to be fitted to the same data ACDM 
already fits so nicely. The hope of a substantial increment of model likelihood was however 
frustrated, so that no real new insight into DE physics has yet been gained in this way. 
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Figure 1: Transfer functions (left) and angular anisotropy spectrum (right) in cosmologies 
with/without coupling and with/without 2 massive :^'s (00/CM models). Coupling and mass are 
selected so to yield an approximate balance. The transfer functions are multiplied by fc^-^, to help 
the reader to distinguish different cases. In the lower frame of the C/ plot we also give the spec- 
tral differences between 00- and CM-models, hardly visible in the upper frame. Here dotted lines 
represent the cosmic variance interval. 



A somehow alternative pattern is an unbiased fit of the scale dependence of DE state 
parameter w^a) to data. But, although restricting to a 2-parameter expression 

w{a) = Wo + {1 — a)w' , (1-1) 

current data hardly do more than fixing a likelihood ellipse provocatorily centered on 
Wo = —1, w' = 0, and up to now, also this approach has failed to give any new insight into 
the DE nature. 

In this paper we however keep within the latter approach, but modifying the parameter 
budget. We aim to test whether existing data are already more constraining, once the range 
of models explored is different. The option we shall explore is suggested by previous results 
based on a Fisher Matrix (EM) analysis [^. Here the authors started with observing 
that spectral distortions due to a coupling between cold dark matter (CDM) and DE 
(parametrized by /3, see below) or to neutrino {v) masses are essentially opposite. This 
holds for both C; and P{k), the CMB anisotropy spectrum and the matter fluctuation 
spectrum (see Figure || for an example). Their EM analysis then assumes that a model 
with vanishing u masses and no coupling (the 00-model) has top likelihood, and tests how 
far one can go, increasing u masses and coupling, keeping within l-cr or 2-a from the 00- 
model. When a CDM-DE coupling is simultaneously allowed, they find that this allows v 
masses well above the limits from the WMAP team , as quoted in Table |l[ 

In this work we shall therefore perform a full Monte Carlo Markov Chain (MCMC) 
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w = —1 


w 




WMAP5 


< 1.3 eV 


< 


1.5 eV 


WMAP5+BA0+SN 


< 0.61 eV 


< 


0.66 eV 



Table 1: Summary of the 2-a (95% C.L.) constraints on the sum of ly masses from WMAP 
5-year, Baryonic Acoustic Oscillations and SuperNova data sets. 



analysis of available data, simultaneously opening the degrees of freedom of v mass and 
CDM-DE coupling. 

Let us remind that CDM-DE coupling jl5U eases the coincidence problem. Within 
ACDM, at z ~ 10^ the DE-matter density ratio is ~ 10"^. With coupling, instead, 
it could already exceed 10^^. With (quasi-) vanishing masses, however, such coupled 



models clash against data, setting a limit (3 < 0.075 [|l5|, [ig, |17 |. 

In turn, at least one u mass eigenstate or, possibly, two of them exceed ~ 0.055 eV (di- 
rect or inverse hierarchy). This follows from solar |^] and reactor Q neutrino experiments, 
yielding Amf 2 — 8 x 10~^eV^, and atmospheric |^ and accelerator beam Q experiments 
yielding Am^ 3 ~ 3 x 10~^eV^. However, the neutrino oscillation experiments do not 
provide us with any information on the absolute scale of masses. 



Cosmology is sensitive to u masses. Already in 1984 Valdarnini & Bonometto |10| 
derived the transfer function for mixed DM models, with DM partially made of massive z^'s. 
Mixed models were widely tested in the Nineties, i/'s then filled the apparently unescapable 
gap between ^lom and ■ That gap is now neatly defined and filled by DE, as already 
outlined. But z^'s becoming non-relativistic cause so strong spectral distorsions that even 
a small contribution to the density budget from them can be tested (for a thorough review 
on effects of massive z/'s on cosmological observables, see ||ll|). This gives the cosmological 
limits on the absolute scale of v masses in Table 1, an order of magnitude more stringent 
than limits from tritium /3-decay experiments. 

The ambitious aim of this paper is then to show that, as both limits can be substantially 
relaxed, (mildly-) mixed coupled models could really be an alternative to the minimal 
ACDM cosmology with the advantage of easing both fine tuning and coincidence paradox. 

The plan of the paper is as follows: In Section |2| we outline the expressions and 
potentials related to the interacting CDM~DE fiuid. In Section ^ we review coupled DE 
(cDE). The data sets and statistically methods used are presented in Section ^. In Section 
m we present our results. Finally we summarize our findings and conclude in Section ^. 

2. Self— interaction potentials 

We shall assume that DE is a scalar field (p, self-interacting through the potentials 

F((/.) = A"+V(/>" RP 

or 

1/(0) = (A"+V0") exp(47r (^V^^p) SUGRA 
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Figure 2: State parameter and its variation in uncoupled RP models. The plot is for h = 0.7, 
rit = 0.046, nc = 0.209. 

admitting tracker solutions. Uncoupled RP yields a slowly varying w{a) state param- 
eter, steadily below -0.85 for A/GeV < 10~^ (see Figure §). On the contrary, uncoupled 
SUGRA yields a fastly varying 'w{a), even faster than any expression ( |1.1D , as is shown 
in Figure ^. Coupling is however an essential feature that we shall be considering and, in 
the next Section, we shall see how it modifies these behaviors. 

Independently of the presence of coupling, for any choice of A and a these potentials 
yield a precise DE density parameter Q^e- Here we use A and Q^e as free parameters; the 
related a value is then suitably fixed. 

Both RP and SUGRA potentials were initially introduced to ease the fine tuning 
problem. In fact, the present DE density can be tuned on today's CMB temperature Tq, 
reading po^de ~ (10 Tq)^; tracker solutions then require a present DE field (po ~ rup so 
that an approximation l^(</'o) ~ Po,de yields mp(lOTo)^ ~ A^"*"". For reasonable a's, A is 
therefore of the order of the energy scale of SUSY breaking or EW transition. 

Fits to data however show that uncoupled RP fails to meet CMB anisotropy data, 
unless the energy scale A stands below ~ 10~^GeV [0] (previous approximate relations 
indicate then a ~ 0.5); this is coherent with the fact that w{a) is then steadily below -0.85 
(~ 95% C.L.) as established in Q for constant w{a) . 

In the case of uncoupled SUGRA, however, while the best fit still leads to low A's, a 
value A ~ 0(0.1 GeV) is within ~ 2(t's from the best-fit model. Let alone the SUSY and 
EW scale, this is still an energy close to the confinement scale. Coupling however causes 
quantitative modifications to these behaviors, that we shall see in the next Sections. 

Quite in general, however, the very fact that the natural representation to describe 
DE is not the one where particle numbers are diagonal, is related to the smallness of the 
mass of the quanta. Such low particle masses, although being a natural consequence of the 
potential expressions, are considered a hidden fine-tuning by some researchers. 
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Figure 3: State parameter in models with 'w{a) given by eq. (1.1) and in a uncoupled SUGRA 
model with A = 0.1 GeV. Wo and w' values (in the frame) selected to yield a behavior close to 
SUGRA, by requiring similar high-z plateau and w{0). Although renouncing to a full coincidence 
at z = 0, the fast variablility of w{a) in uncoupled SUGRA cannot be met by any polinomial w{a). 
The plot is for h = 0.7, fib = 0.046, = 0.209. 



3. CDM-DE coupling 

In the cDE scenario, as for dynamical DE, the scalar field (j) yields a cosmic component 
unsuitable to clustering and showing a negative pressure in the present epoch. Its energy 
density and pressure however read 

p = Pk + ¥{(!)), p = pk-V{^), with = 0V2a^ ; (3-1) 

dots indicate differentiation in respect to r (conformal time), the background metrics being 

ds^ = a^{T) [dr^ - d\^] with d\^ = dr^ + r^{de^ + cos^O dcj)^) . (3.2) 

These expressions show that two regimes are possible. If pk ^ V, the DE state parameter 
approaches +1 {stijf matter) so that DE energy density rapidy dilutes during expansion 
(p oc a^^). In the opposite case V ^ pk, the state parameter approaches -1 and DE is 
suitable to explain the observed cosmic acceleration. 

Let us now consider the possibility that DE is coupled to other components. Inter- 
actions with baryons, constrained by observational limits on violations of the equivalence 
principle (see, e.g. [18|), are almost fully excluded. Similar constraints, however, do not 



exist for DE-CDM interactions. The only constraints then derive from cosmological data. 

The simplest possible coupling is a linear one. It can be formally obtained by per- 
forming a conformal transformation of Brans-Dicke theory (see, e.g., [^]), where gravity 
is modified by adding a (j)R term to the GR Lagrangian {R: Ricci scalar). Coupling causes 
an energy transfer between CDM and DE, so allowing DE to have a non-negligible density 
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even when its state parameter is ~ +1. However, being then dominant, the transferee! 
energy is soon diluted. 

A so-cahed (/>-matter dominated period then occurs when, because of the power leaking 
towards DE, CDM density declines more rapidly than a~^. The increase of <f) then brings it 
to approach and V{(j)) to exceed p^. DE dilution then stops and DE eventually exceeds 
DM density. 

The overall picture is however quite natural. All tenable cosmological models do require 
a dark sector, split into two components with different state equations. The fact that their 
interactions with baryonic matter is just gravitational, leads to requiring that 

^(c)^M^^^r(de)M^^O (3.3) 

(here T^u'^^^ are the stress-energy tensors of CDM and DE, let their traces read T^'^''^'^^), 
while the assumption that CDM and DE are two separate fluids leads to take C = in the 
relations 

T^<i^) z;.^ = + CT^^U^^ (3.4) 
H.^^ = - CT^'U,. , (3.5) 

describing the most general form of linear coupling (incidentally, this shows why DE cannot 
be linearly coupled to any component with vanishing stress-energy tensor trace). 

Assuming two separate fluids is clearly an extra assumption and, if we allow for C 7^ 0, 
when the metric is (|3.2D, these equations yield 

4> + 2-<j) + a^Vl = + Ca'^pc (3.6) 
a ^ 

Pc + 3-pc = - Cpc(t) (3.7) 
a 

Pc being CDM energy density. 

General covariance requires C to be a constant or to evolve as a function of (f) itself. 
Here, instead of C, we shall mostly use the adimensional parameter 

/3 = (3/167r)^/2^pC7 (3.8) 

and consider constant (3 values 0(0. 1) (corresponding to C ~ l/2mp) which, as we shall 
see, meet observational data. 

The above dynamics naturally causes a w{a) behavior significantly different from the 
uncoupled case. In Figure ^ we plot the z dependence of the DE state parameter for the 
SUGRA potential. For the sake of comparison the uncoupled behavior is also shown. 

In the Figure, the transition from a stiff-matter behavior to w ~ — 1 is clearly shown. 
Such transition is later for greater coupling. At low 2, for low /?'s, the uncoupled behavior is 
reapproached and, for j3 as low as 0.05 , the dependence of w on z, when approaching z = 0, 
is even steeper than in the uncoupled case. For greater couplings, however, the behavior 
gradually softens and, for (5 ~ 0.1, the effective value of w at low z is systematically closer 
to ACDM than in the uncoupled case. 
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Figure 4: State parameter in coupled SUGRA model with A = 0.1 GeV, h = 0.7, = 0.046, 
ric = 0.209, as in Figure ^. The uncoupled behavior given there is reported also here. 

Figure § also shows that the scale dependence of the state parameter exhibits a peculiar 
feature, even in the presence of mild couplings {e.g., for /? = 0.05), shifting from negative 
to positive values in a potentially observable redshift range. This behavior is not restricted 
to SUGRA potentials; it holds for RP and other potentials as well. Even approximate 
measurements of the scale dependence of the state parameter, extending up to z ~ 10-30, 
should therefore be able to set direct limits to the coupling intensity. 

The behaviors are qualitatively similar in the case of a RP potential. Figure ^, similar 
to Figure ^, includes a coupling /? = 0.1. It shows, first of all, that w still keeps fairly 
constant values. Furthermore, at low z, the state parameter w is not necessarily increasing 
with z, as it may occur also for the SUGRA potential and is shown in Figure ^ (for 
A = 10~^ GeV), e.g. when (3 = 0.05. As for SUGRA, at low z, w{a) is mostly smaller in 
the presence of coupling. The w = —0.85 limit is bypassed here around A = 10~^ GeV. 
We therefore expect that a fit of coupled RP with data will be already fair for A values 
greater and physically more significant. 

Before concluding this Section, it may be worth defining the coupling function f{(f)), 
through the relation 

C„) . ,3.) 
so that CDM energy density scales according to 

Pc{a) = PoAao/affi^) . (3.10) 
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Figure 5: State parameter in coupled RP for various A values, in the case f3 = 0.1 . Here h = 0.7, 
rif, ~ 0.046, Qc ~ 0.209, as in Figure ^. The uncoupled behavior at z = given there is reported 
also here. 



Then, if we set V = V + pc , the cp eq. of motion takes the (standard) form, 

a ^ 



(3.11) 



as though CDM and DE were decoupled, once the effective potential V is used. 

The CDM evolution (3.10) is then faster than in the absence of coupling. In turn, 
the effective behavior, obtainable by using the potential V, mimics a phantom-like state 
equation, yielding a DE density increase with a, as we would find for w < —1 . 

This makes reasonable that neutrino mass limits are relaxed in a cDE context, as 
they are in the presence of phantom DE. This option, however, does not lead to requiring 
unconventional physics. On the contrary, by assuming quite a general behavior within the 
dark sector and adding the ingredient of u mass, if we are naturally led to high /3 values, 
this will also ease the coincidence problem. 



4. Data and methods 

In this Section we present data sets and methods applied in our analysis. 

To constrain the CMB power spectra, we use the 5 year data from the WMAP satel- 
lite (WMAP5). We use both temperature and polarization data, and calculate the 
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likelihood of our models using the Fortran 90 code provided by the WMAP team^ . 

For the matter power spectrum, we use the results from the 2dF galaxy redshift survey 



21[] . Constraints on the recent expansion history of the Universe is also given by the SNla 
observations from the Supernova Legacy Survey (SNLS) . 

In some cases we also apply additional priors on the Hubble parameter and the baryon 
density of the Universe. The prior on the Hubble parameter of h = 0.72 it 0.08 is taken 
from the Hubble Space Telescope (HST) Key Project IQ. Analysis of the Big Bang 
Nucleosynthesis gives a prior on the baryon content of ^bh^ = 0.022 ± 0.002. 

Using many data sets will of course give stronger parameter constraints, but every 
additional data set might also introduce new systematics. Accordingly, we considered three 
different combinations of data sets in our analysis; (i) WMAP5 data only, (ii) WMAP + 
2dF, (iii) WMAP + 2dF + SNLS + HST and BBN priors (henceforth named "all data"). 

For the MCMC analysis we use the publicly available code CosmoMC ||25| which, 
in turn, uses theoretical power spectra computed by a modified version of CAMB |26|. 
CMB lensing effects are included, to ensure accurate results when comparing with CMB 
polarization data. In the MCMC analysis we use the following set of basic parameters: 
{ cjfo, Wc, ^, T, ris. In lO-'^'^yls, A, /?, Mj^ }. Here: uji,,c are the physical baryon and cold dark 
matter density parameters, tob^c = ^b,ch'^, where h is the dimensionless Hubble parameter; 
9 is the ratio of the sound horizon to the angular diameter distance at recombination; r is 
the optical depth to reionization; Us is the scalar spectral index; As denotes the amplitude 
of the scalar fluctuations at a scale of A: = 0.002Mpc~^. In addition we include the sum of 
v masses, M^, = Tirriy, assuming 3 equal mass v's. Let us outline that A denotes the energy 
scale in DE potentials, while (3 is the coupling parameter between CDM and DE. 

All parameters in the braces are given flat priors, unless otherwise is stated explicitly. 
In our MCMC analysis we also marginalize over the Sunyaev Zel'dovich amplitude, as is 
also done by the WMAP team in their analysis. We assume that the Universe is spatially 
flat. The chains are run on the Titan cluster at the University of Oslo. 

In addition to MCMC runs with the full parameter set presented above, we have also 
repeated the analysis without My or without A and /?, to be able to compare the effects 
of the different extensions of the parameter space. 



5. Results 

In this Section we shall discuss the results of the above MCMC runs. In the cases already 
considered in the literature, we reobtain standard results and likelihood C values. In Table 
|2| we report best fit values and 1-a errors for different models. 

When including the /?, A and degrees of freedom we see that all the other parame- 
ters stay within \a shifts from their previous mean values. As could be expected, the error 
bars increase on some of the parameters, especially iVc- 

The most intriguing part of our outputs however concerns My and (3. Already from 
Table I one appreciates that the upper limit on u masses has doubled. This can also be seen 

^http://lambda. gsfc.nasa.gov; version v3 
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Parameter 



ACDM + z/'s 



w = const. cRP + i/'s cSUGRA + z/'s 



WMAP only all data all data all data all data 



2.244 2.258 2.247 

± 0.066 ± 0.061 ± 0.062 



2.260 2.260 
± 0.061 ± 0.065 



0.1156 0.1098 0.1132 0.1039 0.1042 

± 0.0078 ± 0.0040 ± 0.0069 ± 0.0062 ± 0.0084 



10^^ 



1.0401 1.0401 1.0402 1.0401 1.0406 

± 0.0030 ± 0.0030 ± 0.0030 ± 0.0029 ± 0.0030 



0.085 0.087 0.085 0.087 0.088 

± 0.017 ± 0.017 ± 0.017 ± 0.016 ± 0.017 



(eV) 
(95% C.L.) 



< 1.20 



< 0.66 < 0.94 



< 1.13 



< 1.17 



(95% C.L.) 



<0.17 



<0.18 



logio(A/GeV) 
(95% C.L.) 



< -4.2 



< 6.3 



0.955 0.962 0.958 0.969 0.970 

± 0.017 ± 0.014 ± 0.015 ± 0.015 ± 0.018 



ln(10^M3) 



3.053 3.045 3.049 3.055 3.057 

± 0.043 ± 0.040 ± 0.040 ± 0.040 ± 0.041 



0.691 0.713 0.711 

± 0.075 ± 0.056 ± 0.059 



0.723 0.717 
± 0.062 ± 0.069 



Ho (km/s/Mpc) 



67.0 
± 4.4 



70.1 
± 2.1 



69.7 
± 2.2 



71.8 
± 2.5 



71.9 
± 2.7 



-2 ln{C) 1329.39 1407.25 1407.38 1407.44 1407.33 

Table 2: Best fit values and 1-a error bars. In all fits we allow for u masses. The first 9 lines concern 
primary parameters. Only upper limits on M^, (3 and A are shown. These variables are discussed 
more thoroughly in forthcoming 2-D plots. Likelihood values are almost model independent. 
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in the Figures ^ and where we show marginahzed and average hkehhood distributions 
on the most specific parameters of this work: (3, log(A/GeV) and M^. 

Figures from the top to the bottom hne refer to fits based on increasingly wide data 
samples. Taking WMAP5 data only, the marginalized likelihood distribution gives just a 
95% C.L. upper hmit (3 < 0.28 and (3 < 0.23 for SUGRA and RP potentials, respectively. 
When low-z data are simultaneously considered, the marginalized likelihood distributions 
exhibit significant maxima, which could be naively interpreted as a ~ 2 a detection of 
f5. A physical reason of the effect can be found in the actual tension between ojc mean 
values, obtained from pure WMAP data or including low-z data, visible also in Table |2|. 
A significant DE amount at high z could reduce there the required value of ujc by 1-2%, 
and set it within 1 a from its all-data value. 

In turn, models with /? are somehow more likely, although one must be cautious on 
this point, where non-Gaussian behaviors become important, in the presence however of 
tiny signals. It is true, in fact, that the peak of the marginalized likelihood, in some cases, 
exceeds 0-0 models by almost 2 it's. But the peak is not so pronounced among average 
likelihood values: here the top value is atmost double, in respect to 0-0 models. Finally, 
if we consider the top overall likelihood, it does not exceed ACDM likelihood. This seems 
to conflict with the fact that the likelihood of 0-0 models in RP and SUGRA cosmologies 
equals the likelihood of w = —1, among w = const cosmologies. But, of course, the 
likelihood of individual models can easily behave differently from averaged likelihoods. 

As far as the A scale is concerned, the 2-a upper limit on log^^Q A is much softened 
compared to studies of uncoupled models. The usual upper limit ~ 0.1 GeV, for SUGRA 
models, shifts now above lO^GeV. Something similar occurs for the scale A in RP models, 
which is now consistent with data, within ~ 2(t's, up to ~ 10~^GeV. This had been 
somehow predicted from Figure ^, as outlined there. 

Finally, in the M,^ distributions, the gray line shows the distributions in the absence 
of coupling, and allows to appreciate why the 95 % upper limits in Table |2| have almost 
doubled. The maxima in the likelihood distributions are far less accentuated here, than 
for p. They are somehow stronger in the SUGRA than in the RP case. It is however clear 
that there is no hint of j/-mass detection in these plots. 

In all these plots there are discrepancies between average and marginalized likelihood 
distributions. They are particularly relevant as far as P is concerned. 

Such discrepancies, first of all, are a safe indication of non-Gaussian distributions. In 
the case of /? they can be better undestood in association with some Figure herebelow, and 
they are surely the basic reason to cast serious doubts on the formal /3 detection from the 
marginalized likelihoods. The non-Gaussian behavior is minimal for masses. 

The most significant plots of this paper are Figures |^ and |^, concerning the SUGRA 
and RP potentials, respectively. In these Figures, curves yield 1- and 2-a contours for 
marginalized likelihood distributions. On the contrary, colors yield average likelihood dis- 
tributions. The plots concern the f3-M,j and the /3-logiQ(A/GeV) planes. 

These plots allow, first of all, a better undestanding of the discrepancy between mean 
and marginalized likelihood distributions. Let us consider the correlations between (3 and 
log]^o(^/G'eV), shown in the lower panels. They indicate that the same (3 values are consis- 
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Figure 6: Marginalized (solid black line) and average (dotted black line) likelihood of cosmological 
parameters in SUGRA models. For M^we also show the corresponding likelihood distributions 
obtained in the case of a standard ACDM-|-M^model (gray lines). 
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Figure 7 : Marginalized (solid line) and average (dotted line) likelihood of cosmological parameters 
in RP models, taking into account the whole set of available data. 
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Figure 8: Two parameter contours for the SUGRA model. Solid lines are 1- and 2-cr limits for 
marginalized likelihood. Colors refer to average likelihood, and the 50% likelihood contour from 
the average likelihood is indicated by the dotted line. 

tent with a fair range of A's. When marginalizing over the other parameters, A included, 
this increases the weigth of high /3 values. 

In particular, Figure |8| exhibits a peculiar high-A tail: these 2~D plots show that, for 
suitable couplings f3 ~ 0.1 , models with A up to 10^^ GeV (close to GUT) are allowed 
within 2-(t's. No such features is present for the RP potential. Again, however, one 
appreciates that A ~ MeV is now allowed for RP if /5 ~ 0.15-0.17 . 

Thus, when it comes to the fine-tuning paradox, the SUGRA potential keeps a more 
satisfying solution than the RP potential, as a A scale close to the EW scale or to the scale 
of the soft SUSY break is consistent with data. Also RP, however, does no longer yield 
just unacceptably low energy scales. 

One should however also consider these potentials independently of the supposedly un- 
derlying physics, as examples of rapidly or slowly varying w{a). Using such potentials we 
could actually inspect the behaviors within these extremes, making recourse to a single pa- 
rameter A, instead of using, e.g., the parametrization Wq-w', which requires 2 independent 
variables. 

Let us then point out that, in the RP case, closer to w = const., both the marginalized 
and the average likelihood distributions on (3 exhibit a maximum. 

Also in the left panel of Figure |^ the likelihood peak is quite far away from /? = 0, 
Mjy = 0. In the marginalized likelihood distribution, such model is within the 95% C.L.. 

More in general, from the left panels of Figures ^ and showing the 2D likelihood 
distributions in the [3-Mi, plane, we want to stress the following points: 

(i) The solid curves in the SUGRA case substantially overlap FM results from in the 
SUGRA case. 

(ii) Marginalized and average likelihood are different, as is expected in the presence of a 
non-Gaussian behavior, but not significantly so. 

(iii) Top likelihood models are found for non-vanishing (3 and M^; the best fit values are 
13 ~ 0.07 and ~ 0.35 eV for both potentials. 
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Figure 9: As previous Figure, for a RP potential. 

(iv) Models with P up to 0.2 and Mi^ up to 1.4 are allowed with both potentials within the 
95% C.L. . Again a result apparently potential independent. 

(v) The limits on /3 and are strongly correlated, as expected from the FM analysis Q. 

(vi) According to the marginalized likelihood distributions, (3 = Mi, = models are allowed 
for both potentials. In the RP case, however, they lay outside the 1-a limits. 

(vii) The fact that upper limits on v masses are loosened by almost a factor 2 arises from 
the degeneracy between f3 and M^. 

This confirms the results from |5| and shows that one should be cautious when ex- 
tracting neutrino mass limits from cosmological observations, as they can heavily depend 
on the assumed model range, as is the case here. 



6. Summary and conclusions 

One of the major discoveries of physical cosmology is the existence of the dark sector of 
the Universe. Particles belonging to the standard model of elementary interactions account 
for not more than 5 % of its energy budget. Observations then go farther and show that 
the remaining 95 %, which interacts with standard model particles just through gravity, is 
not a single component, but needs to be modeled at least through two independent fluids, 
with drastically different state equations. 

All existing data can be accomodated in a scheme assuming no energy exchange be- 
tween these two dark components. This is simple, but leads to the well known fine tuning 
and coincidence paradoxes. 

In this work we have explored the option that, in the dark sector, energy exchanges 
occur and may be described through a linear and constant coupling. Clearly, this is just the 
next approximation beyond assuming zero exchange and is however a phenomenological 
approach. Hopefully, it may lead to constraints helping to single out a precise theory, just 



as is done when the expression (LI) is taken for the state parameter scale dependence. 

Such option was already explored in the past. The limits found for the CDM-DE 
coupling were then rather deceiving. For /3's large enough, the coincidence paradox could 
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Log a 



Figure 10: Evolution of the density parameters in a SUGRA model with coupling and v mass. 
Colors refer to different components, as specified in the frame; lines to the different models: M^—^, 
(3 = (continuous line); M„=0.5 eV, (3 = 0.085 (dotted hne); M„=l.l eV, (3 = 0.17 (short dashed 
line); M^=1.2 eV, /3 = 0.22 (long dashed line). 



be significantly attenuated; but the allowed /? range did not allow to go so far in this 
direction. 

The critical ingredient of this work amounts to considering simultaneously coupling 
and u masses. Massive z^'s would be a further component of the dark side, but unlike 
from CDM and DE, however, they are particles already known from the standard model. 
Furthermore, their role appears unessential to the formation of cosmic structures, although 
their mass can actually modify the matter distribution over the largest scales. 

When one allows simultaneously i' mass and f3, we find models with significantly higher 
likelihood. The upper limits on both parameters are softened and increase by an essential 
factor ~ 2. This is a critical factor in both fields. The new allowed M^, values approach 
the i^-mass detection area in forthcoming tritium decay experiments like KATRIN p7|| . 
Moreover, when considering an external prior on M,^ from earth-based experiments, the 
strong degeneracy between the coupling parameter /3 and the neutrino mass Mi, can be 
broken, gaining new insight into the DE nature p^ . 

Simultaneously, the new allowed P values open the possibility of a critically modified 
DE behavior. In Figure 10 we show the scale dependence of the cosmic components for 
various f3~My pairs. They tell us that models are allowed, within the 95 % C.L., for which 
DE is still ~ 1 % of the cosmic energy budget up to z ~ 10^; at higher redshifts it decreases 
just because the photon-neutrino fraction increases, while we approach matter-radiation 



equality. Such fraction attains 2-3 % for the most coupled case we considered in Fig. 10 
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corresponding to C = 0.9/mp with Mu=1.2 eV; this is at 2.3 a's from the best fit. Let us 
remind that ACDM cosmologies prescribe that, at 2; ~ 10^, DE bears less than 1 : 10^ of 
the total energy density. 

Available data do not yet force us to require a non-vanishing CDM-DE coupling. The 
statistical analysis of data still leads to an intricate situation, where marginalized and 
average likelihoods exhibit discrepancies. Furthermore, the likelihood distributions on the 
coupling (3 exhibits some dependence on the selected self-interaction potential. Using a RP 
potential, the 0-0 option appears rather unlikely, both through marginalized and average 
likelihood distributions. A SUGRA potential, instead, yields a higher likelihood for the 0-0 
option. Accordingly, we believe that current data do not allow a claim of "/3-detection" , 
while they certainly allow to put upper limits to the coupling, which can be so large, within 
the 95 % C.L., to yield C = l/2mp. 

Furthermore, it is fascinating to notice that, if shortly forthcoming particle data will set 
a lower limit to Mi^, in the range they are allowed to explore, this will imply an almost model 
independent CDM-DE coupling detection, opening the way to a deeper understanding of 
the dark sector of the Universe. 
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